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Abstract — We develop an approximate analytical technique 
for evaluating the performance of multi-hop networks based 
on beacon-less CSMA/CA as standardised in IEEE 802.15.4, 
a popular standard for wireless sensor networks. The network 
comprises sensor nodes, which generate measurement packets, 
and relay nodes which only forward packets. We consider a 
detailed stochastic process at each node, and analyse this process 
taking into account the interaction with neighbouring nodes via 
certain unknown variables (e.g., channel sensing rates, collision 
probabilities, etc.). By coupling these analyses of the various 
nodes, we obtain fixed point equations that can be solved 
numerically to obtain the unknown variables, thereby yielding 
approximations of time average performance measures, such 
as packet discard probabilities and average queueing delays. 
Different analyses arise for networks with no hidden nodes and 
networks with hidden nodes. We apply this approach to the 
performance analysis of tree networks rooted at a data sink. 
Finally, we provide a validation of our analysis technique against 
simulations. 

I. Introduction 

IEEE 802.15.4 [1], a popular standard for the physical 
layer and medium access control for low-power wireless 
networks, underlies the now well-known ZigBee standards for 
wireless sensor networking. The standard defines two types of 
CSMA/CA algorithms - slotted and unslotted, and is named on 
the basis of the algorithm chosen - beaconed and beacon-less, 
respectively. With the growing importance of wireless sensor 
networks in industrial applications Q, we need analysis and 
design techniques for multi-hop IEEE 802.15.4 networks. In 
this paper we develop a new approximate analytical technique 
for multi-hop beacon-less IEEE 802.15.4 networks. 

In Figure [T] we depict the queueing schematic of a network 
of the type that we are concerned with in this paper. There 
are sensor nodes, namely, 1, 3, 4, 5, 7, and 8, that generate 
measurements (at rate A;) that need to be transported to the 
base station (BS). The sensor nodes can also serve as relays 
for the traffic of other sensor nodes, e.g., Node 3 is a relay for 
nodes 4 and 5. Two additional nodes, 2 and 6, serve only as 
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A 4 A 5 



40 As Os 




Fig. 1 . A WSN consisting of nodes arranged in a TREE topology. Source 
node i has packet generation rate A^. Relay nodes 2 and 6 forward the data 
of their predecessors to the base-station BS. 



relay nodes, due to the limited range of the radios associated 
with the sensor nodes. The nodes use unslotted CSMA/CA, 
as standardised in IEEE 802.15.4, to contend for the wireless 
medium and transmit their packets. Even the stability analysis 
of such networks is a difficult problem, due to the fact that 
the service rates applied to queues are non-monotone with the 
queue occupancy of the other queues. In principle, the entire 
network can be modeled via a large coupled Markov chain, 
with the state at each node being the number of packets in its 
queue, and the state of the contention process for the head-of- 
the-line (HOL) packet. Such an approach is well known to be 
intractable even for a network with a single contention domain 
and saturated queues (see Q, Thus all researchers have 
taken recourse to developing approximate analyses. Table [I] 
summarizes some of the recent work on modeling approaches 
for CSMA/CA networks where we have listed the major 
limitations of the proposed models as well. Owing to the 
need for low power operation and large coverage, in general, 
a sensor network is multi-hop with the presence of hidden 
nodes. Further, for networks that carry measurement traffic, the 
queue occupancy varies with time. In addition, unsynchronized 
operation is preferred to save large overheads involved in the 
synchonizations of devices. Unlike other models, we consider 
all of these attributes in our modeling approach. 



TABLE I 

A COMPARISON OF LITERATURE ON ANALYTICAL PERFORMANCE MODELS OF CSMA/CA NETWORKS. 



Authors 


Network Scenario & Limitations 


Modeling Approach 


Kim et al.(5] 


beacon-less 802.15.4; rare packet ar- 
rivals, star network; no hidden nodes; 
no queueing at nodes; with ACKs 


2-dim discrete time Markov chain for each node with states: (backoff stage, backoff 
counter) while in backoff, and separate states when transmitting and idling; with 
unknown transition probabilities; coupling the per node models to get the unknowns 
via fixed point approach 


Singh et al.(6J 


beaconed 802.15.4; (un)saturated, star 
network; no hidden nodes; with ACKs; 
infinite queue size; homogeneous nodes 


Markov renewal process at each node with cycle length dependent on the number of 
nodes available to attempt in the first backoff period of the cycle; transition probabilities 
and conditional expectations using node decoupling; coupling via stationary probabilites 
of remaining nodes; saturation analysis results used for finite arrivals 


Qiu et al.|7J 


802.11 DCF; (un)saturated, star net- 
work; with hidden nodes and ACKs; 
infinite queue size 


discrete time Markov chain with states as the set of transmitting nodes having unknown 
transition probabilities, together with RSS measurements between each pair of node to 
obtain CCA and packet failure probabilities; iterative procedure for solving unknowns 


He et al.|8J 


beaconed 802.15.4; saturated, star net- 
work; no hidden nodes; with ACKs; 
homogeneous nodes 


two Markov chains for each node: one embedded at the end of transmissions in the 
common channel, and the other during backoffs embedded at each slot; both having 
states as (backoff stage, backoff counter); chains coupled to get a fixed point approach 


Martalo et al.|9J 


beaconed 802.15.4; unsaturated, star 
network; no hidden nodes; no ACKs; 
finite queue size 


two semi-Markov processes: one for tagged node having unknown transition probabili- 
ties, and the other for the shared radio channel that couples all the per node models to 
obtain the unknown quantities; MGFs used for finding queueing delays 


Goyal et al.HOI 


beacon-less 802.15.4; unsaturated star 
network; no hidden nodes; no queueing 
at nodes; with ACKs 


tracked the number of nodes with non-empty queues assuming that the probability of m 
nodes having non-empty queues at any given time is same as the probability of m — 1 
empty nodes getting a new packet to send while a non-empty node is sending its current 
packet (for coupling individual node models); assumed that the CCA attempt processes 
of other nodes with non-empty queues are Poisson processes (decoupling) 


Lauwens et al.| 1 1 1 


beacon-less 802.15.4; saturated, star 
network; no hidden nodes; no ACKs 


semi-Markov model for each node with unknown and non-homogeneous transition 
probabilities based on the backoff stage of the node; also found the distribution of 
backoff intervals and iterated over these together with the stationary probabilities of the 
semi-Markov processes 


DiMarco et al.1121 


beacon-less 802.15.4; unsaturated, 
multi-hop network; no queueing at 
nodes; with hidden nodes and ACKs 


3-dim discrete time Markov chain for each node having state: (backoff stage, backoff 
counter, retransmission counter) with CCA and collision probabilities (unknown) as 
state transition probabilities (decoupling); assumed independence of node processes to 
find expressions for the CCA and collision probabilities using stationary probabilities 
of different Markov chains (coupling) 


Shyam-Kumarl 13 1 


beacon-less 802.15.4; (un)saturated, 
multi-hop network; with hidden 
nodes and ACKs; infinite queue size; 
mismatching results 


saturated analysis using continuous time Markov chain with states as the set of 
transmitting nodes having unknown transition probabilities; found CCA and failure 
probabilities using steady state probability of the Markov chain; unsaturated network 
modeled as a Markov renewal process constituted by the nodes having non-empty 
queues, saturation analysis results used in each cycle with the set of non-empty nodes 



Our Contribution 

We consider a multi-hop WSN consisting of sources and 
relays arranged in a tree topology operating in the beacon-less 
mode with ACKs. Each node has an infinite buffer space and 
may operate in the saturated or unsaturated regime with fixed 
packet length. Different analysis techniques are developed for 
networks with hidden nodes and networks devoid of hidden 
nodes due to the difference in activity lengths perceived by a 
node in the presence and absence of hidden nodes. We model 
the stochastic process evolving at a node by incorporating the 
influence of the other node in the network by their (unknown) 
time averaged statistics. After making certain approximations, 
we use renewal theoretic arguments to couple individual node 
processes and obtain the unknown variables as fixed points of 
an iterative scheme. Although this decoupling (or mean-field) 
approximation is popular in such situations, our more detailed 
model incorporating several issues not considered together 
hitherto requires a careful handling of the analysis. We identify 
and calculate two QoS measures for each source node, viz., 
the packet delivery probability and end-to-end packet delay, 
in terms of these variables. We observe that in a multi-hop 
network, the packet delivery probability falls sharply before 
the end-to-end packet delay becomes substantial. 

The rest of the paper is organized as follows: In Section [H] 



we give an overview of the unslotted CSMA/CA mechanism. 



Section III explains the node behaviour and elaborate upon 
the different analysis techniques used for networks without 
hidden nodes and networks with hidden nodes. Section [IV] 
compares our analysis with simulations. We conclude the 
paper in Section [V] 

II. Unslotted CSMA/CA for Beacon-less IEEE 
802.15.4 Networks 

Node behaviour under unslotted CSMA/CA is shown in 
Figure [2] A node with an empty queue remains idle until it 
generates a packet or receives one from its predecessor nodes. 
When a node has some data to send (non-empty queue), it 
initiates a random backoff. The first backoff period is sampled 
uniformly from to 2 macmmBE — 1, where macminBE is 
a parameter fixed by the standard. The backoff period is 
specified in terms of slots where a slot equals 20 symbol times 
(T s ) and a symbol time equals 16 /is. The node then performs 
a CCA (clear channel assessment) to determine whether the 
channel is idle. If the CCA succeeds, the node does a Rx- 
to-Tx turnaround, which is worth 12 symbol times and starts 
transmitting on the channel. The failure of the CCA starts a 
new backoff process with the backoff exponent raised by one, 
i.e., to macminBE+1, provided it is lesser than the maximum 
backoff value given by macmaxBE. The maximum number 



of successive CCA failures for the same packet is governed 
by macMaxCSMABackoffs, exceeding which the packet is 
discarded at the MAC layer. The standard allows the inclusion 
of acknowledgements (ACKs) which are sent by the intended 
receivers on a successful packet reception. Once the packet is 
received, the receiver performs a Rx-to-Tx turnaround, which 
is again 12 symbol times, and sends a 22 symbol fixed size 
ACK packet. A successful transmission is followed by an 
InterFrame Spacing{lFS) before sending another packet. 

When a transmitted packet collides or is corrupted by the 
PHY layer noise, the ACK packet is not generated which is 
interpreted by the transmitter as a failure delivery. The node 
retransmits the same packet for a maximum of aMaxFrameRe- 
tries times before discarding it at the MAC layer. After 
transmitting a packet, the node turns to Rx-mode and waits for 
the ACK. The mac AckWaitDu ration determines the maximum 
amount of time a node must wait for in order to receive 
the ACK before concluding that the packet (or the ACK) 
has collided. The default values of macminBE, macmaxBE, 
macMaxCSMABackoffs and aMaxFrameRetries are 3,5,4 and 
3, respectively. 

III. Analytical Model 

As shown in Figure [T] the network consists of sensors, 
relays and a base-station (BS). We consider time invariant 
links and static (immobile) nodes. We assume that the sensors 
generate traffic according to independent point processes. 
These constitute the aggregate external arrival process for 
the sensor network. Each sensor transmits its data to the 
next-hop node according to the topology (Tree or Star). The 
intermediate nodes along a route may be relays in which 
case they simply forward the incoming traffic, or they may 
be sensors which transmit their own packets as well as the 
received packets. Based on the network congestion, the nodes 
may discard packets due to consecutive failed CCAs or frame 
retries. Figure [2] shows a typical sample path of the process 
evolution at a node i. 



The exact analysis of such a process is intractable. We, 
therefore, make several simplifications: 

(51) All packets are of the same fixed length. Hence, the 
DATA transmission duration is fixed. 

(52) The internode propagation delays are of the order of 
nanoseconds, and hence can be taken to be zero. 

(53) Zero packet capture. As a result, all the packets involved 
in a collision get corrupted. 

(54) The IFS is neglected. It means that the receiver does not 
discard the packet received during the IFS. 

(55) The packet error probability at each link is fixed and 
known. As already mentioned, the links are not time varying 
and we can empirically obtain the error characteristics of a 
link. 

(56) ACK packets are short and, therefore, not corrupted by 
PHY layer noise. 

(57) A node's CCA succeeds when there is no transmission 
by any node in its CS range at the time of initiation of its CCA. 
Recall that in the standard, the channel state is averaged over 
the 8 symbol duration. 

(58) The time taken by a transmitting node (called the trans- 
mission period and denoted by T) for the activities of success- 
ful transmission and collision are the same. If the transmitted 
data collides at the receiver, the macAckWaitDuration for the 
transmitting node is equal to the sum of the turnaround time 
of 12 T s and the ACK duration of 22 T s , a total o/34 T s (see 
Figure 

(59) All nodes have the same CS and transmission ranges. 
This implies that if a node i is in CS (resp. transmission) 
range of another node j, then so is node j in the CS ( resp. 
transmission ) range of node i. 
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Fig. 2. Process evolving at Node i where it goes through various states. 

Since the process evolutions at various nodes are coupled 
and each node (source or relay) can be present in any state, the 
random process representing the system is multi-dimensional 
and non-homogeneous. Further, different nodes can interact 
with different subsets of nodes based on their Carrier Sense 
(CS) range and positions, the process is also non-symmetric. 



Fig. 3. The transmission period T which includes the DATA transmission 
time T x , the turnaround time of 12 T s , and the MAC ACK duration of 22 T s 
(or the macAckWaitDuration of 34 T s , by (S8)). 

Further, we make use of a decoupling approximation 
whereby we model each node separately, incorporating the 
influence of the other nodes in the network by their average 
statistics and as if these nodes were independent of the tagged 
node. 

The channel activity perceived by any node i is only due 
to nodes that are within its CS range (denoted by fii). Let 
us consider the durations during which the queue at node i 
is non-empty. During these times, node i alternates between 
performing CCAs, and transmitting the packets from its queue. 
Let the CCA attempt process at node i be a Poisson process 
of rate fa conditioned on being in backoff periods. Let (3j 
be the CCA attempt rate of node j in 0,^ over all time, and 
approximate the CCA attempt processes for the nodes in f2j 
by independent Poisson processes. By modeling simplification 



(SI), we assume that all packets entering node i have the same 
fixed length, and hence take the same amount of time when 
being transmitted over the medium (denoted by T). As a result, 
we observe that as perceived by node i in conditional time, 
its transmission periods and the intervening backoff periods 
constitute a renewal process, as shown in Figure [4] Note that 
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Fig. 4. The renewal process obtained by observing the process at a node i 
after removing all the idle time periods from the original process shown in 
Figure 2 The renewal epochs are denoted by {-X,- } and the cycle lengths 
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Fig. 6. Node j 6 Qi finishes its backoff, performs a CCA, finds the 
channel idle and starts transmitting the DATA packet. Node i finishes its 
backoff anywhere in the shown 12 T s duration and there is no other ongoing 
transmission in Qi, its CCA succeeds and it enters the transmission duration. 
As a result, the DATA packets may collide at r(i) £ Qj and/or r(j) S 
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Fig. 5. (a) expands a non-empty period (conditional time) and displays a 
renewal cycle when there are no hidden nodes in the system, and (b) shows 
the expanded non-empty period when hidden nodes are present. 

f3j is equal to the product fijbjqj where fjj is the CCA attempt 
rate of node j in backoff periods, bj is the fraction of time 
node j is in backoff provided its queue is non-empty, and qj 
is the fraction of time node j's queue is non-empty. Figure [5] 
shows a renewal cycle of node i. The dashed transmission 
durations belong to the nodes in fi^. 

Before discussing the difference in activity lengths in the 
two network types (as shown in Figure [5J, we list the various 
packet collision scenarios. Firsty, the majority of packet colli- 
sions in a wireless network occur due to the presence of hidden 
nodes. However, there are scenarios where transmissions from 
two transmitters placed within the CS range of each other 
can overlap, known as Simultaneous Channel Sensing and 
Vulnerable Window and are explained in Figures [6] and [7] 
respectively. Note that we are making use of (S2) to avoid 
including insignificant timing details here. 

Further, the ACK packet size is just a fraction of DATA 
packet size (e.g. compare an ACK packet (22 symbols) with a 
DATA packet of length 260 symbols at PHY layer). It means 
that if a node is in its transmission period T, then we can 
say with high probability that it would be in the T x (DATA) 
duration. This fact further says that the probability of a node's 
CCA initiating in the vulnerable window of another node is 



Fig. 7. Vulnerable window of 4 T s in the transmission period of node j £ Q; 
during which node i's CCA attempt would be successful. 

negligible. Hence, we can safely assume that the CCAs seldom 
fail due to ACKs. Similarly, the small ACK packets can rarely 
cause packet collisions. 

We now return to elucidate the difference in length of 
activity periods in a renewal cycle. 

(a) Absence of Hidden Nodes: When a node transmits, all 
the other nodes in the network can hear it resulting in CCA 
failures for other nodes that try to assess the channel in 
the transmission period. Also if two nodes are involved in 
simultaneous channel sensing, the activity period may extend 
from T to a maximum of (T + 12 T s ). Since 12 T s « T, 
we can assume that the activity period is only a single 
transmission period of duration T. 

(b) Presence of Hidden Nodes: Since a node's hearing ca- 
pacity is limited, it may not perceive the activities of all the 
nodes in the network which may cause dilation of activity 
period. As shown in Figure [H] nodes j and k are hidden from 




Fig. 8. Dilation of transmission period as perceived by node i. 

each other with respect to the receiver node i. Let node j 
starts transmitting while nodes i and k are silent. Since node 
k cannot hear node j, it can start its own transmission while 
the transmission from node j is still going on. For node i, 
this looks like a dilated transmission period, whose length is 
denoted by 7T (e//) . 

Since each node has a different set of neighbours, we 
make the following approximation to simplify our analysis: 



The mean length of Dilated Activity Period as perceived by 
node i is equal to the length of mean busy period of an 
M/D/oo queue where the deterministic service time is equal 
to a single transmission period T, and the arrival process is 
approximated as a Poisson process having rate equal to the 
aggregate transmission initiation rate of all the nodes in the CS 
range of node i. Note that the above assumption is equivalent 
to saying that all nodes in fij are hidden from each other, and 
hence results in a larger mean for the dilated activity period. 

Let Tj be the rate of successful CCAs by node j, or 
equivalently, the rate of transmission initiations by node j, 
which is equal to the product j3j(l — ctj) where ctj is the 
probability of CCA failure for node j. For any node i, let the 
aggregate transmission initiation rate for nodes in fij is 



(i) 



Assuming node j £ fi, has started transmission and another 
node k £ n ( starts at time u < T, the expression for T< e "> 
can be written recursively as 



T, 



( e ^) = Te-^ T + f^(u + T} em )Cie-^du 



A. Topologies with No Hidden Nodes 

By the assumptions made so far, we have a superposition 
of independent Poisson processes during the backoff periods 
of node i. Define the following: 



Pi 
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and 



9i = 



= 1 - e 
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Since ZigBee is designed to operate at fairly low packet 
generation rates at which packet discard probability (due to 
successive CCAs or frame retries) is very low, all the recursive 
equations written below do not take packet discards into 
account. 

1) Probability of CCA Failure, af. The probability of CCA 
failure is the probability of the occurrence of at least one 
transmitting node in the CS range of node i, given that node 
i performs a CCA. Defining N t (t) and N^ F \t) as the total 
number of CCAs and number of failed CCAs in the interval 
(0,t], respectively, we have: 



ai = lim — „ = lim 



jv| f) (t) 



Applying Renewal-Reward Theorem (RRT) lfl4ll here, we get 



Nj 
IV", 



(3) 



where is the mean number of failed CCAs and Ni is 

the mean number of total CCAs in a cycle. The mean renewal 
cycle length, Wi, is given by 



Wi 



9i + xhT + (1 - m) Cl T + (1 - th)(1 - Ci)(T + W l ) 



By taking the reward to be the mean number of total CCAs 
by node i in a renewal cycle, i.e., JV», and applying RRT, we 
can write 

Ni = th + {l-r)i)ci + (l- 77,0(1 -Ci)(ftT + JVi) 

Note that node i performs CCAs at the rate of fa for the 
entire transmission period T of another node, i.e., a total of 
faT CCAs all of which fail and the cycle continues after this 
period. Again, let the reward be the total number of failed 
CCA attempts by node i in a renewal cycle, say, N£ , which 
can be written using (S7) as 

NP = (l-^Xl-c^ftT + iVf)) 

On rearranging, the CCA failure probability is given by 

r) i + (l-r)i)c i + (l-'rii)(l-Ci)PiT 

2) Packet Failure Probability, ji: A transmitted packet can 
fail to be decoded by its intended receiver due to a collision, 
or due to noise. We do not take packet capture into account 
here (see (S3)). The probability of packet failure, 7^ is defined 

as: 



7i 



Pi + (1 - pi)k 



(5) 



where is the probability of packet collision and is the 

probability of data packet error (known by (S5)) on the link 

between i and its receiver r(i) due to noise. Recall that ACKs 

are not corrupted by PHY noise (see (S6)). Defining MAt) 
(c) 

and Ml (t) as the total number of transmissions and total 
number of collisions in interval (0, t], we can write 



7» 



Urn 

t— ^00 



Mf\t) 
Mi{t) 



1 



lim 

t— >oo 



Mf\t) 
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Dividing each term by t and applying RRT here as before, we 
get 



M 



(C) 



M 
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Mj 
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1 - 



Wi 
Mj 
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where Ml is the mean number of collided packets and M, 
is the mean number of transmitted packets in a cycle. If we 
ignore packet discards, there is a single transmission by node 
i in each renewal cycle, Mi becomes unity. We denote the set 
of nodes that can cause interference in successful reception at 
r(i) by I r u), which is composed of two sets: 



C 



r(i) - {3 G ^ '■ 3 S n * and 3 G / r( 4 )} 



r (2) __ 



{j £,yV :j i^li and j £ I r{l) } 



such that 



r (2) _ 



Note that the receiver r(i) is in C r X (see Figure 




9). 



Fig. 9. Interference region around the receiver. R( cs ) is the CS range which 
is assumed to be equal for all nodes by (S9). 

In this case, the collisions can happen only due to the 
simultaneous channel sensing event. Let the reward be the 

mean number of collisions of node i's packets in a renewal 

(c) 

cycle, say, M( , which can be written as 



M 



(C) 



m (l - cxp{-12 T s J2 f3 3 }) + (i - my 

r(C) 



+ (l- m ){l- Ci )Ml 

Note that the first two terms are simultaneous channel sensing 
events. On rearranging ,the packet collision probability at node 
i is given by 



n 



(l - cxp{-12 t s J2 Pj}) + (i - Vi)d 



1-(1-^)(1- Cl ) (6) 

The packet failure probability can then be calculated using 
Equations [5] and [6] 

B. Topologies with Hidden Nodes 

For each node i, we now need to keep track of channel 
activity in its CS range only. Further, even if a node j £ 17^ 
finishes its backoff first when there is no activity in f2 i; its 
CCA need not be successful because it may be blocked by 
another node in the system hidden to i. Hence, we redefine rji 
and gi as 

ft . 1 



ft + E 



and g l 



ft 



The probability of CCA failure, a,, can be evaluated on 
similar lines as before with the exception that T gets replaced 
by T { ; n \ i.e., 



(l-TJiXl-C;)/^ 



(e//) 



r/ 4 + (l-^)((l-c l )/3 l T I (e//) +c ! ) 



(7) 



2) Packet Failure Probability, 7$: We first need to find 
the probability that nodes in a given set are not transmitting 
using a product approximation due to the unavailability of 
joint distribution of processes. Let hi be the fraction of time 
a node is not transmitting while its queue is non-empty. The 
reward Hi in this case is the mean amount of time node i is 
not transmitting in a renewal cycle given by 

H g. + q-^Xl-cQTf^ 

l-(l-7Ji)(l-Ci) 

which in turn provides 

h . = * = gi + (l- Vl )(l- Cl )Tt ff) 

W > gi + mT+{l-m)mT+{l- m ){l-Ci)T^ ff) 

The fraction of time node i is not transmitting (unconditional) 
is denoted by hi and is equal to 



hi = (l-qi)+qihi 



(8) 



The expression for collision probability can be found as 
follows: 

1. The first term accounts for the fact that Node i started 

transmitting in the presence of at least one transmission by 

(2) 

the hidden nodes set, i.e., C , \. 

r(i) 



R 



(i) 



m ( 1 - Yl h i 



ipr <2) 



2. The second term accounts for the scenario when Node i 
started its transmission as a simultaneous channel sensing with 
a node in f2j in the presence of at least one transmission by 



the hidden nodes set, i.e., C, 



(2) 

w 



R 



The expression for Cj remains the same. 

1) Probability of CCA Failure, a*: The mean renewal cycle 
length Wi can be written as 

Wi = gi + mT + (1 - m ) Cl T + (1 - - Cl )(Tj eff) + Wi) 

Note that we have replaced T in the last term by 
because this period is susceptible to dilation. 



f = (l-ttte/i- n h A 



3. The third term considers the case where Node i starts 
transmitting in the absence of any ongoing transmission by 
a hidden node but it encounters a simultaneous transmission 
by a node in C^i) anywhere in the corresponding 12 T s period 

(2) 

and/or a transmission by a node in C\!y i.e., a hidden node 
anywhere in its activity period T. 



Rl ' = Vi I rL ec (2)^ \f j [ 1-exp ^ -12 T s . ^ ' j 

SpP (1> 




E 



exp < — T 



4. The final term says that Node i started its transmission as 
a simultaneous channel sensing with a node in C^X in the 
absence of any hidden node. This event surely ends up in a 
collision at the receiver r(i). 



E 



R 



(4) 



n ^ 



All the other possible scenarios have a negligible probability. 
The overall packet collision probability is given by 



r)i + (1 - J7 4 )ci 



(9) 



The packet failure probability, ji, can now be calculated using 
Equations [9] and [5] 

C. Evaluation of Other Quantities 

1) Average Service Rate (erj: Each packet that reaches 
the HOL position in the queue at a node can have multiple 
backoffs and transmissions before it is successfully transmitted 
or discarded. We define Zi as the average time spent in 
backoff, and Yi as the average transmission time until the 
packet is successfully transmitted or discarded at the MAC 
layer. Then the average service rate, cr, at node i can be 
calculated as: 



1 



Z, + Yi 



(10) 



Using the default values from the standard, the mean backoff 
durations can be calculated (shown in Table [II]). Define the 
following quantities: 

B, = (70+20+158a i +318a 2 +318af+-318a*-12o^) 

T,(l) _ I (70 + 20X1- 



(i-°?) 



I 248014(1-0^) . 566o|(l- Qi ) 

' T, — r- 



(l-c?) 



(1- 



884a-?(l- Qi ) 1202o:|(l- 
+ T, ,5, + 



tP> = (78+158+318+-318+318) 



where Bi refers to the mean backoff duration until the packet it 
transmitted or discarded due to successive CCA failures, T^ l) 
has the interpretation of the mean backoff duration given that 

(2) 

the packet transmission was successful, and T> ' is the mean 
time spent in backoff given that the packet was discarded due 



Backoff 


Mean Backoff Duration 


Mean Backoff Duration 


Exponent 


with successful CCA(T 3 ) 


with failed CCA(T S ) 


3 


70 + 20 


70 + 8 


4 


150 + 20 


150 + 8 


5 


310 + 20 


310 + 8 



TABLE II 

Mean backoff durations in symbol times, T b 



to successive CCA failures. Then the quantities Z t and Yi can 
be calculated as: 

Zi = a=T i (2) +(l-a^)[T/ 1) +7 1 (QfT/ 2) +(l-Q 5 ) 

[T^ 1) +7i(Q^T I (2) +(l-a=)[T i (1) +-7 1 (a=Tf ) + (l-Q=)T i <1) )])])] 
Yi = (l-af)[T+ 7i (l-a?)[T+7i(l-c«f)[T+ 7i (l-a i i )T]]] 

2) Aggregate Arrival Rate (vi), Goodput {Bi) and Discard 
Probability (Si) for a Node: The arrival process at each node 
consists of packets which are generated at the same node (if 
it is a sensor) and the packets to be forwarded. The aggregate 
arrival rate at a node i can be written as: 



(11) 



keVi 



where \ is the packet generation rate at sensor node i, Vi 
is the set of all its in-neighbours, and 9i is the rate of packet 
transmission by the node which are successfully received at 
r(i) (known as goodput). An enqueued packet at a node is 
successfully received by the intended receiver in a manner 



shown in Figure 10 
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Fig. 10. The goodput of a node is denned as the rate of successfully received 
packets by its receiver. 

Assuming that the queueing system is stable and has a 
steady-state solution, 9{ is given by 



^(1 - Si) 



(12) 



where Si is the probability of discarding a packet due to 
consecutive CCA failures or successive failed retransmission 
attempts, and is given by 



(13) 



Note that if the queue at node i is saturated, then the goodput 
9i is equal to <7j. Note that while calculating <jj, we have 
considered the packet discards into account. 

3) The node non-empty probability, qii To find the ex- 
pression for Qi, assuming that all the arriving packets reach 
the HOL position (i.e., no tail drops) and applying Little's 
Theorem, we get 



<7, 



(14) 



Further, for a saturated node, the quantity ^ is equal to 1. 

4) Obtaining bi and fy: In order to find bi, the fraction of 
time a node is in backoff provided it is non-empty, we embed a 
renewal process in conditional time where the renewal epochs 
are those instants at which the node enters the random backoff 



period after a packet transmission or packet discard. We use 
the RRT to find the expression for 6, as 



bi = 



Bi 



Bi + (1 



)T 



(15) 



The quantity fa, the rate of CCA attempts in backoff times, 
is based on the backoff completion times irrespective of the 
transmission attempt of the packet since a packet retransmis- 
sion is considered the same as a new packet transmission. 
Thus, for the calculation of fa, it suffices to observe the 



process at node i only in the backoff times. Figure 11 shows 
the residual backoff process where after completion of the 
backoff duration, the node performs a CCA. Using the result 




CCA mkyl'cJs and node docs a R\-n>-T\ uaa:an>ui:d 



Consecuiivc CCA i'adnrcs and nodi' discards llie packet 



Fig. 11. Evolution of backoff periods in conditional (backoff) time. The y- 
axis is in units of a backoff slot (20 T s ). The range indicates that the random 
backoff duration is uniformly distributed within it. 

of Kumar et al. [4 | with their collision probability replaced by 
our CCA failure probability, on, the expression for fa can be 
written as 



fa = 



1 



Bi 



(16) 



D. End-to-End Delay Calculation 

We use the two parameter Whitt's QNA |[T5l to calculate 
the mean sojourn time at each node based on an approximate 
Moment Generating Function (MGF) of service time. We can 
write an expression for service time in a recursive manner by 
allowing infinite CCAs and retransmissions. To account for 
the discarded packets, we allow only a fraction of transmitted 
packets to join the next-hop neighbour's queue based on the 
discard probability at that node. We denote the service time at 
node i by Si and let Bi denotes the length of random backoff 
duration, which is assumed to have an exponential distribution 
with rate fa for node i. Then, Si can be written as 



Bi 
Bi 
Bi 



Si 
T 
T 



S, 



w.p. 
w.p. 
w.p. 



(1 
(1 



7*) 



where Si is a random variable with the same distribution as 
Si. The MGF of B h denoted by M Bz {z), is equal to 
Therefore, we can express the MGF Mg t (z) of service time 

Si as 



fa(l - ai )(l - ji)e' 



z + fa{l - ai)(l - 7ie-* T )' 

The first two moments of the service time, E(5j) and E(Sf ) 
can be calculated by differentiating the MGF M^fz). The 



QNA procedure commences from the leaf nodes and converge 
to the base-station in a sequential manner. We first calculate 
the squared coefficient of variance of service time, denoted by 
Cg., at each node, i.e., 

E(^ 2 ) 1 
5 * (E(S 4 )) 2 

Ignoring the packet discards at the MAC layer, the net arrival 
rate Aj at a node i is given by: 

!Aj + Aft for source-cum-relay nodes 
ker< 
y Afc for relay nodes 

For any node i, let pi be defined as 



Pi 



AjE(Si) for leaf nodes 
A,*E(>£j) for internal nodes 



We now calculate the squared coefficient of variance for the 
departure process {c 2 D .) for a node which uses the squared 
coefficient of variance for the interarrival times at that node 
(c A .). Here we include the probability of a packet discard 
which is equal to <5j. 

c\ = i + *i(^(4 4 -i) + (i-p?)(^-i)) 

Further, c\. is calculated as 



4 



Finally, the mean sojourn time at a node i is given by 

PiE(S i )(c 2 Ai +c 2 Si 



A, 



E(Si) (17) 



2(1 -ft) 

The end-to-end mean packet delay for a source node j, 
provided that the set of nodes along the path from this node 
to the BS is Lj, is given by 



(18) 



E. Packet Delivery Probability for each Source Node 

The fraction of packets generated at each source node 
that reaches the BS is called the packet delivery probability, 
denoted by pf for each source node i. Let the set of nodes 
constituting the path from a source node i to the BS be Li. 
Then, assuming that the drop events are independent from 
node to node, the expression for pf el ^ turns out to be 



P., 



(del) 



(19) 



where Sj is the packet discard probability at node j. 

F. Global Iteration Scheme 

The global iteration scheme is shown in Figure [12] We start 
with a vector {0, 0, Xi} corresponding to {a,, Vi\ for each 
node i, and repeat the procedure until the quantities converge. 



Fig. 12. The global iteration scheme. The square boxes indicate the derived 
quantities and the round boxes indicate the known quantities. 

IV. Numerical and Simulation Results 

For the verification of our analytical model, we use QualNet 
(v4.5) simulator with the default parameter values and a fixed 
packet size of 70 bytes. However, the QualNet implementation 
is devoid of ACKs. Therefore, we compare results only for 
the ACK-less scenarios (although the analysis is for a general 
scenario). The network shown in Figu re [T3] has no hidden 
nodes whereas the one shown in Figure |16| has hidden nodes 
(the dependency graph is also shown). Note that all the nodes 
are sources (some of which also serve as relays) with identical 
packet generation rates that are simultaneously increased. 

In networks with no hidden nodes, we have accurate results 
for various quantities, viz., CCA failure probability (aj, 
packet discard probability (S), packet failure probability (72), 
goodput (8) and node non-empty probability (q). Note that 
the discard probability is plotted on a log-log scale. It can 
be observed that the queue at node 1 gets saturated when 
the packet generation rate reaches 20 pkts/sec. However, the 
goodput f?i starts to drop drastically after 5 pkts/sec. Moreover, 
since node 5 has a lower aggregare arrival rate than that of 
node 1 (v 5 < v{), its queue remains unsaturated even when 
qi becomes unity. The goodput #5 also shows departure from 
the input rate ^5 at around 5 pkts/sec. 

In the network with hidden nodes, the results for 6, 7, 9 
and for q match well. It can be seen that the goodput starts to 
deviate from the input rate at around 2.5 pkts/sec for node 1 
and around 4 pkts/sec for node 8. There is a mismatch in CCA 
failure probability (a) which is reflected in the plots for b, (3, 
and A. However, even for this scenario, the analysis captures 
the overall qualitative behavior. If the network is operated only 
until the drop probability is less than 10%, then the errors in 
all the measures for each node are less than 17%. 

Further, as observed for both the network types, the node 
packet discard probability increases to an impractical value 
before the average node delay becomes substantial, e.g., node 
1 discards 2-5 packets for every 1000 packets when the 
average node delay is 5-6 msec. Hence, the overall delivery 
probabilities for the sources act as performance bottlenecks 
for these networks. 



Fig. 13. A 10 nodes Tree topology. Note that each node can hear every other 
node. 
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Fig. 14. Results for Node 1 for the tree topology of Figure \\3\ 
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Fig. 15. Results for Node 5 for the tree topology of Figure [T3] 
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Fig. 16. A 10 nodes TREE topology. Also shown is the dependency graph 
where the dotted lines connecting two nodes indicate that the nodes are in 
CS range of each other. 
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Fig. 17. Results for Node 1 for the tree topology of Figure [T6| 
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Fig. 18. Results for Node 8 for the tree topology of Figure [T6] 



V. Conclusion 

We have developed an approximate stochastic model for 
the performance analysis of beacon-less IEEE 802.15.4 mul- 
tihop wireless networks, with arrivals. The model permits 
the estimation of several average performance measures. Our 
model is accurate for topologies devoid of hidden nodes and 
works reasonably well for topologies with hidden nodes in 
terms of the packet discard probability, failure probability and 
throughput. We calculated the mean end-to-end delays and 
packet delivery probabilities for each source in the network, 
which gives an idea of the QoS provided by the sensor network 
at a specific packet generation rate at source nodes. The results 
suggest that, for the relatively small size tree networks that we 
have studied, to operate in the low-discard low-delay region, 
the packet arrival rates at nodes should not be greater than a 
packet every few seconds (e.g. a packet inter-generation time 
of 5 to 10 seconds at the sources). 
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